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Abstract. We examine the entropy of stationary noncquilibrium measures 
of boundary driven symmetric simple exclusion processes. In contrast with 
the Gibbs-Shannon entropy [Tl 1101 . the entropy of noncquilibrium stationary 
states differs from the entropy of local equilibrium states. 



1. Introduction 

In the last decade important progress has been accomphshed in the understand- 
ing of nonequilibrium stationary states through the study of stochastic lattice gases 
([HIH] and references therein). 

The simplest nontrivial example of such dynamics is the one-dimensional simple 
symmetric exclusion process on the finite lattice {1,2, ...,A^ — 1} with particle 
reservoirs coupled to the sites 1 and — 1. In this model the microscopic states are 
described by the vector 77 = (77(1), ?7(2), . . . , r]{N — 1)), where = 1 if the site i 
is occupied and 77(1) = if the site is empty. Each particle, independently from the 
others, perform a nearest-neighbor symmetric random walk with the convention 
that each time a particle attempts to jump to a site already occupied the jump is 
suppressed. At the boundaries, particles are created and destroyed in order for the 
density to be a at the left boundary and /3 at the right boundary, < a, /3 < 1. 

We denote by //^^ the stationary state of this system which is a probability 
measure in the space of configurations and which can be expressed in terms of a 
product of matrices [5]. Since the particle number is the only conserved quantity 
in the bulk, in the scaling limit N — >■ 00, i/N x € [0, 1], the system is described 
by a single density field p{x), x G (0,1). The typical density profile p{x) is the 
stationary solution of a partial differential equation with boundary conditions. In 
the context of symmetric exclusion processes, 

p{x) = a{\ — x) -\- j3x . 

The nonequilibrium stationary states exhibit long range correlations [13j which 
are responsible in the large deviations regime for the non locality of the free energy 
functional O [2] ■ More precisely, if 7 stands for a density profile different from the 
typical one p, the asymptotic probability of 7 is exponentially small and given by 

where the so called nonequilibrium free energy Va^p is a non local functional. 
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Since in equilibrium the probability of such large deviations is determined by the 
induced change in the entropy, it is natural to investigate the entropy of nonequi- 
librium stationary states. 

Denote by &n{i^^) the Gibbs-Shannon entropy of a state ly^: 

e^(^^) = -^/^^Wiogi^^(77), 

■n 

where the sum is carried over all lattice configurations rj. Recently, Bahadoran [T] 
proved that for a large class of stochastic lattice gases the Gibbs-Shannon entropy of 
nonequilibrium stationary states has the same asymptotic behavior as the Gibbs- 
Shannon entropy of local equilibrium states. In our context of boundary driven 
symmetric simple exclusion processes this result can be stated as follows. Denote 
by g the product measure 



N-l 



1=1 

Thus, at site i, independently from the other sites, we place a particle with prob- 
ability p{i/N) and leave the site empty with probability 1 — p(i/N). Bahadoran 
proved that 

The long range correlations of the nonequilibrium stationary state is therefore not 
captured by the Gibbs-Shannon entropy. 

Derrida, Lebowitz and Speer jlOj showed that for the symmetric simple exclusion 
process the difference 

converges as N ^ oo, and that the limit depends on the two points correlation 
functions. Hence, the long range correlations appear in the first order correction to 
the Gibbs-Shannon entropy. 

In this article we examine the entropy of the stationary nonequilibrium states 
fJ,^ In the classical Boltzmann-Gibbs theory of equilibrium statistical mechanics 
[12], the steady state p^ (?/) of a niicrostate 77 is given by 

f^piv)^^^exp{-m{v)) (1.1) 

where /? is the inverse of the temperature, H{ri) the energy of 77 and Zn(/3) the 
partition function. The Boltzmann entropy is then defined as the limit, when the 
degrees of freedom N of the system converges to infinity, of 1 /N times the logarithm 
of the number of microstates with a prescribed energy: 



S{E) = lim lim log ( 1{\H {j]) - N E\ < SN}) , 

V 

where the summation is performed over all configurations rj and where 1{j4} is the 
indicator of the set A. The pressure P{P) is defined by 

P(/3)== hm -i-logZ^(/3) 

JV-s-oo iV p 
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and the Boltzmann entropy is related to the pressure function by 

SiE) = inf {/3P(/3) + I3E} . 

I3>Q 

In view of (jl.ip and by analogy, we define the energy of a microstate 77 as 
— log /i^^ (77) and the entropy of the stationary nonequilibrium measure A*^/3 by 

S^AE) = lim lim N-'\og( ^ 1{\N-' log ^,^ArJ) + E\ < 6}) . 

We propose in ()2.5p a variational formula for the entropy function Sa,p in terms 
of the nonequilibrium free energy Va^p and the equilibrium Gibbs-Shannon entropy, 
that we conjecture to be valid for a large class of boundary driven stochastic lat- 
tice gases. This formula is based on a strong form of local equilibrium, stated as 
assumption (H). We present in (|2.9p an explicit formula for the entropy function 
Sa,p and wc show in ()2.13p that it is strictly concave, being the Legendre transform 
of a strictly concave function Pa,p, identified as the nonequilibrium pressure. This 
last point is proved in sectional 

In Section [3] wc compute the entropy of stationary nonequilibrium measures of 
boundary driven zero range processes and in Section [7] we show that the entropy of 
the nonequilibrium stationary states is different from the entropy of the local 
equilibrium states p- In Section O wc determine the energy band and describe 
the density profiles with lowest and largest energy. In Section |6l we examine the 
isentropic surfaces and in the appendix we show that the strong form of local 
equilibrium holds for the symmetric simple exclusion process by using the ideas of 
i- 



2. Stationary nonequilibrium entropy function 

Fix an integer iV>l, 0<q;</3<1 and let An ■= {1, . . . , - 1}. Denote 
by := {0, 1}"^" the configuration space and by 77 the elements of Hn, so that 
ri{x) = 1, resp. 0, if site x is occupied, resp. empty, for the configuration rj. We 
denote by a^ yj] the configuration obtained from 77 by exchanging the occupation 
variables ri{x) and rj{y), i.e. 

{■q{y) li z = x 
•nix) if z = y 
ri{z) if z ^ x,y, 

and by a^r] the configuration obtained from 77 by flipping the configuration at x, 
i.e. 



{a^v){z) 



1 — i]{x) if z ~ X 
r]{z) if z ^ X. 



The one-dimensional boundary driven symmetric exclusion process is the Markov 
process on fl^ whose generator Ln can be decomposed as 

Ln — Lo,N + L-^N + L+^N , 
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where the generators Lq.n, L-^n, L+,n act on functions / : Q.n — >■ M as 

N-2 

x=l 

{L^,Nf Kv) = {a[l - r,(l)] + (1 - a)7il)} [/(a^r,) - /(r,)] 

(L+,Ar/)(r,) = {/3[l-ry(7V-l)] + (l-/3)rK^-l)}[/(a^-S)-/(77)] . 

We denote by r/t the Markov process on fijv with generator L n . Since the Markov 
process rjt is irreducible, for each A'^ > 1, and < a < /3 < 1 there exists a unique 
stationary state denoted by /i^^. 

The entropy function Sa,i3 '■ 1R+ {^cxd} U [0,log2] associated to the nonequi- 
librium stationary state A*^^ is defined by 

liii- , , 

5^0 N^oo N 



S^,p{E) = hmjhn^ilog 1 { |iV-i log + i?] < <5} 



whenever the limits exist. To keep notation simple, we sometimes denote Sa^p by 
S. 

Note that we may include in the sum IJ,^ p{Ti): 

S^AE) = E + Km lim ^\og ^i%{\N-Hog ^Ji%{r^) + E\ < 5} . (2.1) 
In particular, 

Jo^AE) = E - So^AE) (2.2) 
is the large deviations rate function of the random variables —N~^ log /i^^ {rj) under 
the probability measure /x^^. 

At equilibrium a = /?, the stationary state /i^ ^ is a Bernoulli product measure 
with density a and the entropy function is given by 

So^iE) := S^AE) = ~"\ ) , (2.3) 

V log a — log(l — a) / 

where 

s{9) = log + (l-0)log(l- 61) 
represents the Gibbs-Shannon entropy. This formula is valid for E in the energy 
band [E^{a), E^{a)] where 

_E_(a) = — log { max(a, 1 — a)} , £'+(«) = — log { min(Q;, 1 — a)} . 

In the case a = 1/2 the energy band is reduced to the point log 2 and 5i/2(log2) = 
log 2. Outside the energy band we have Sa{E) — ~oo. 

Identity (|2.3p can be derived from the large deviations principle for the random 
variable —N~^ log ^{r]), which in the equilibrium case is an average of i.i.d. 
random variables. 

Denote by ( • , • ) the scalar product in L^([0, 1]). Let A4 be the set of measurable 
profiles m : [0, 1] [0, 1] equipped with the topology induced by weak convergence, 
namely m„ — > to in if and only (m„, G) — > (m, G) for every continuous function 
G : [0, 1] — ?► R. For every m ^ A4 the nonequilibrium free energy [9l[2] 14,.^ (m) of 
TO is defined by 

' m(x) 1 — TO(a;) F'{x)' 



/ r ■m[x) L — m[x) f [x)^ 

= I |m(.x)log-p^+(l-TO(..))log^-^ + log^|dx 
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where F e C^([0, 1]) is the unique increasing solution of the non linear boundary 
value problem 

^"-("^-^)F(r^' (2.4) 
_F(0) = a , F{1) = /3 . 

To keep notation simple we frequently denote Va^p by V. 

Decompose the set A^r into r = e^^ adjacent intervals Ki, . . . ,Kg-i of size eN 
and denote by M = {Mi, . . . , Mr) the number of particles in each box. Let 

i^iMi, Mr) = E i{ E ^(^) = ^'^1' ■ • • ' E ^(•^) = ^^'■} 

be the probability to find Mj particles in the interval Kj, 1 < j < r. Denote by 
/j,^^(-|M) the probability measure fi^ ^ conditioned to have Alj particles in Kj, 
j = 1, . . . ,r. The set of configurations rj such that X^xG-R" '?(^') ^ -^^i' 1 ^ J ^ is 
denoted by J7Ar(M) and its cardinality by Zjv(M). We shall assume that for every 
< a < /? < 1, 



1 

lim lim sup sup sup — 

£->-0 Ar_>oo M rjenwCM) N 



log{Z^(M)M^^^(r;|M)} =0. (H) 



We present in the appendix a formal derivation of this hypothesis. 

Assumption (H) states that the stationary state ^ conditioned on the number 
of particles on macroscopic intervals is uniformly close in a logarithmic sense to 
the uniform measure as the number of intervals increases. As we shall see, this 
alternative formulation of local equilibrium plays a central role in the investigation 
of the entropy of stationary nonequilibrium measures. The first main result of this 
article provides a variational formula for the entropy function. We claim that for 
every 0<a</3<l, £'>0, 



where 



Sc,jj{E) = sup {S(m) : Va,^p{m) + §(m) = E} , (2.5) 

S{m) = - f . 

Jo 



s{m{x)) dx 



This formula is a straightforward consequence of assumption (H) and the large 
deviations for the nonequilibrium stationary state A*^^. Indeed, we may rewrite 
l^a ^s l^a i3iv\^)'^(^)- Hence, by definition of the entropy and by assumption 
(H), 

S(E) = lim lim lim —log V lUN-Hogiy(M) - N-Hog Zn{M) + E\ < s\ . 

The previous sum can be rewritten as 

E^iv(M)l||iV^Mogi^(M)-A^-MogZAr(M)+£;| < (5| . 

M 
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Recall that N''^\ogiy(M) - -V{Tn) and that 7V"MogZjv(M) S(m) where m(-) 
is the macroscopic profile associated to M: 



m = 

i=l 



, p, = M,/iNe), = {[Nx^], . . . ,[Nx,+i] - 1} . 



Since for a fixed s the sum over M has only a polynomial number of terms in N and 
since Z^CM) is exponentially large in N, only the term which maximizes Zjv(M) 
contributes. The result follows. 

Note that by [3l Theorem 4.1], the functional Va,i3 + S is continuous in A4. 

2.1. The nonequilibrium pressure. Let A : — > R+, P : R ^ M be given by 

f/^ dx 
Ja [x<^ + [1 - x^]'^' ^2.6) 
P{0) = P^AO) ■■= log {^0^) , MO, 

P(0) = —log 2. As we shall see in p.l3|) . P is the Legendre transform of the 
entropy function S and may thus be identified with the nonequilibrium pressure. 
An elementary computation shows that log A is strictly increasing on the intervals 
(—00,0) and (0,oo) and that lime^±o log ^(0) = ^oo. Moreover, 

_ 1 rWi^-'^))" 1 (l + a:)log(l+a:)-^logx AjO) 
^^"^ ^(^)^2i(a/(i-a))« l + x (xi/Ml)(x«-i+x«)i/e ''^ + ^°g;3_a 

(2.7) 

for 61 7^ 0, and 

/ 1 /■'^ dx \ 
P'{0) = log(- / _ - ) . 

We prove in Section |4] the following properties. 
Lemma 2.1. The junction P is a strictly concave function. Moreover, 
lim {P{e)-eE^] =0, lim {P{e) - OP' {9)) = 0, 



= E+ia,P) = log / — ^- 

nnn|a;, 1 — a;} - 

f I dx \ 

i?_(a,/3) = log(- / — -) 

\ p ~ a Jq, max|a;, 1 — x\ / 



It follows from this lemma that lim6i_>±oo 9{P'{9) — E^i = 0, and that for each 
E G {E-,E+), there exists a unique 9e = 9{a, (3,E) £ R such that 

P'i9E) = E . (2.8) 

Define the functions 7± : R — R+ by 

7+(0) = 7+(«,/3,e) = max{(^)', (y^)'} , 
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and let We = Wa p e ■ [7-(^'b)7 7+(^'-e)] ^ [0, 1] be the monotone function given 

by 

, , I r dt 

We{x) 



Clearly, W£;([q;/(1 — a)]^^) = 0. On the other hand, the change of variables t = 
[x/{l-x)f'' shows that We{[P / {I- P)f'') = 1 in view of the definition of yl(6'). Let 
/iB : [0, 1] ^ [1-{0e),1+{6e)] be the inverse of H^e so that /ii5(0) = - a))"^ , 

hE{l) = {P/{l-P)Y-. 

2.2. An explicit formula for Sa,p- We are now in a position to present an explicit 
formula for the entropy function S. We claim that for every < a < /3 < 1, 

E_{a,f3) <E < E+{a,l3), 



Indeed, consider the variational problem p.Sp . Let 6 be the Lagrange multiplier 
and let R{m, 6) be the function defined by 

R{m,e) ^ S(m) - e{V{m) + S{m)- E} . 

Since by OE] {SV/6m) = log[m/(l-m)]-log[F/(l-i^)], the conditions {5R/5m) = 
dgR = imply that 



|m(a;) log F{x) + [1 - m{x)] log[l - F{x)] - log (|-^) } dx 



(2.10) 



where F is the unique increasing solution of the non linear boundary value problem 
(f2^ . We report the first identity in (|2?T0| to (f2^ to get that 



F" F''-! - (1 - Ff 



F' I (1-F)« + F» 



F' 



■/^^ - fl - z)''-! 

Since dldz \B~-^ log((l - z)'' + z")] = ^ '-^ and (logF')' = ^''/F' we 

(1 — z)'' + 

deduce from the previous equation that 

F' = A[{1-FY + F'']^'' (2.11) 

for some positive constant A determined by the boundary conditions satisfied by 
F: 

F'jx) 

Jo [(l-F)« + F«]i/«"'^- 



The change of variables y = F{x) shows that A = A{9) is given by p.6p . 

Recall the definition of 7±(0). Let gg : [0,1] [7_(0), 7+(0)] be given by 
gg ~ {F/{1 — F))^ and observe that 

g', = A{e)e{gy' + l){gl-'+g<>Y/r 
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Define Ue : [7-(6i), 7+(6l)] ^ [0, 1] by 

and remark that gg ~ Ug^ . 

In the second equation of (|2.10p replacing mhy gg/{l-\- go) and F' by the right 
hand side of identity ()2.1ip . we obtain that 

\f{T^^^>ggg{x)~\og{l+ge[x))]dx - log - . 

Performing the change of variables y = ge {x) , we get that 

1 /■(/^/(i-ffl)'' I a;logx-(l+x)log(l+a:) _ Ajd) 

(2.12) 

In view of the explicit expression for P' , we may rewrite the previous identity as 
P'{e) = E. Therefore, by (g^, 6* = ^b, and hence Ug = We, ge = hs- Moreover, 
in view of (|2.10p . the density profile m which solves the variational problem (|2.5p 
is m = hs/il + hs]- This proves ((2J)) . 

2.3. A variational formula for Sa,i3- We conclude this section showing that 
Pa,p is the Legendre transform of Sa^p and can therefore be identified with the 
nonequilibrium pressure. 

For every 0<Q!<,3<l,i?>0, 

SiE) = M {BE -Pie)} . (2.13) 

If E belongs to the energy band {E^, £'+) the infimum is attained at 6'e given by 
([^ and 

SiE) = esE - Pies) . 

Moreover, SiE±) = and SiE) = -oo ii E ^ [E-,E+\. 

By abuse of notation we shall call S the Legendre transform of P. Usually 
the Legendre transform is defined as a supremum and involves convex functions. 
However, by taking a minus sign we may transform convex functions into concave 
functions and supremums into infimums. 

The proof of (|2.13|) is simple. In section [5] we show that SiE) = — oo outside 
[E_,E+] and that SiE±) = 0. On the other hand, by Lemma [2Jl eE ~ Pie) is a 
monotone non-decreasing function for E > E+. Hence, for E > mie^s.{OE — 
Pie)} = lime^-oo{OE - Pie)}. By Lemma O again, {eE - Pie)} converges to 
—oo, for E > E+, E = E+, respectively. Therefore, by the first observation of 
the proof, SiE) = in{{eE - Pie)} for E>E+. The case E < E^ is analogous. 

By Lemma [2n \img^±oo{eE - Pie)} = +oo for E e iE-,E+) and the function 
e eE - Pie) is strictly convex on R. Hence, mi{eE - Pie)} = es E - Pies), 
where 6'^; solves (j2.8|) . We may rewrite this expression as 0eP'(0_b) — Pi&s)- In 
view of (|2.7|) . to conclude the proof of (|2.13p it remains to show that the first term 
on the right hand side of (|2.7p multiplied by e and computed at 6* = coincides 
with (j2.9|) . This can be shown by performing the change of variables u — hsix) in 
(j2.9p and recalling that We is the inverse of He- 

It follows from ()2.13p that 5* is concave and that P, the Legendre transform of 
the entropy, can be identified with the pressure. 
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The equilibrium case can be recovered by letting a (3. In this case, 

1 £: + log(l-a) 



hE(x) = 



log[a/(l - a)] E + \og{a) 
E + log(l - a) 
E + log(Q;) 



3. Boundary driven zero range processes 

We compute in this section the entropy of stationary nonequilibrium measures 
of boundary driven zero range processes. The model is described by a positive 
integer variable 'q{x) representing the number of particles at site x S Ajv. The 
state space N'^" is denoted f^Ar. At exponential times one particle jumps with rate 
g{r]{x)) to one of the nearest-neighbor sites. The function g : N — > M+ is increasing 
and g{Q) = 0. Wc assume that the system interacts with particle reservoirs at the 
boundary of Ajv whose activity at the right is (p+ > and at the left is > 0. 
The microscopic dynamics is defined by the generator 

Cn — Cq^n + ^-,N + C+,N 

where 

= E [fiT^'^'^'v) - /W] + g^x + 1)) [/(T^+^'^r?) - /(r;)]} 

x=l 

C^,N.fiv) = 5('7(1)) ifiS^v) - fiv)] + V- ifiS+v) - fiv)] , 

c+,Nf{ii) = giviN - 1)) [f{s^_,v) - fiv)] + v+ [f{s^-iv) - fiv)] , 

in which 

{ri{z) if z^x,y 

ri{z) - 1 if z^x 

ri(z) + 1 if z = y 

is the configuration obtained from 77 when a particle jumps from x to y, and 

^ " ' [ 77(2) ±1 if z ~ X 

is the configuration where we added (resp. subtracted) one particle at x. Note 
that, since g{0) = 0, the number of particles cannot become negative. 

The invariant measures of the boundary driven zero range processes can be 
computed explicitly. Let ipN ■ A^v — > be the linear interpolation between <y9_ 
and 

/ X \ X 
(Pn{x) = (^1-— + ^'^+- (3.1) 

The invariant measure is the product measure whose marginals are given 

by 



771 



where Z{(p) = 1 + J2k>i f'^/lgi^) ' ' is the normalization constant. 
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Denote by R : M+ — > K+ the density of particle under the stationary state with 
activity equal to ip on both boundaries: 

and by $ : M+ — > 1R+ the inverse of R: 

$ = i?-^ . 

Under the stationary state, the typical density profile p : [0, 1] — > M+ is the 
unique solution of the elliptic equation 

r A<f>(p) = , 

\p(0) = i?(^_), p(l) = i?(^+), 

where A stands for the Laplacian. As TV t cx), the activity profile tpN introduced 
in p.ip converges to $(p): 

lim ifN = , 

A* — j-oo 

where t^^v : [0, 1] M+ is the function defined by (^Ar(O) = ip-, (Pn{x/N) = (pn{x), 
X G Ajv, <^Af(l) ~ 'P+ and extended to the interval [0, 1] by linear interpolation. 
The weight of a configuration rj under the stationary state m^_ ^ is given by 



N-l 

{ri) = exp^ \r]{x)\ogipN{x) - log[g{l) ■ ■ ■ g{r]{x))] - log Z{lpn{x))^ 



In the special case where g{k) — l{k > 1}, the weight 'm,^_ ^^(?7) is a function of 
the empirical density. Hence, if in this case we define for a profile p : [0, 1] — R-(-, 

n^..vAp) ■■= I p{x)\og^{-p{x)) - log Zmp{x)))dx , (3.2) 



we have that 



In general, the weight of a configuration is not a function of the empirical 
density but a function of the field associated to the variables ^(x) = ri{x) — 
log[g(l) • • • g{'ii{x))]/ log ipN{x). 

The nonequilibrium free energy functional is easy to compute in the context of 
zero range boundary driven systems since the stationary state is a product measure. 
A simple computation shows that 

T/ I ^ [\ u '^(pW) , zmpjx))) ^ 

Jo ^{P[x)) Z{^{p[x))) 

To present an explicit formula for the entropy function in this context, we need 
to introduce some notation borrowed from the theory of large deviations of i.i.d. 
random variables. Fix > 0, let M : M — >■ R.). be given by 



Z{^)^^g{l)-.-g{k) 



where F^(fc) = k- [logip]'^ log[g(l) • ■ ■ g{k)], and let R^{a) = M'^{a)/M^{a). The 
large deviations rate function for the mean of the i.i.d. random 
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variables ^{x) = ri{x) — [logt/j] ^ log[5(l) • • • g{r]{x))], 1 < x < N — 1, distributed 
according to Ti^^^ is given by 

I^ix) = xR-\x)~logM^iR-\x)) . 

In the particular case where g{k) = l{k > 1}, we get that Z{lp) = (1 — 'p)~^, 
R{^) = cp/il-ifi), Hp) = p/il + p), F^{k) = fc, = rKa-), M^{a) = (l-^)/(l- 
ipe"), R^{a) = (/?e"/[l - ^pe"] so that 

X 

iJx) = a; log (1 + a;)log(l + x) - log(l - x) 

<^{x) Z^x)) (3.3) 

X log log , , • 

ip Z[ip) 

We emphasize that formulas p.2p and p.3p have been deduced only in the case 
g{k) = l{fc > 1}, and may not hold in general. 

For each < <y9_ < if+, define the entropy function Sip_,ip^ : M — s- M by 



whenever the limits exist. Wc may introduce in the sum ^^{i]) to get that the 
entropy function is equal to 

E + lim lim ^logm!^ | iV-Mogm^ + ^ ^ 4 

Since log (77) may be expressed in terms of the variables {£,{x) '■ 1 < x < 

— 1}, which are independent under the large deviations principle gives 

that ^ 

S,fi_,^+{E) - E = - M /$(p(j.))(A(x)) , 
^ Jo 

where the infimum is carried over all profiles A : [0, 1] M such that 

nl 

X{x)\og<S>{p{x)) -\ogZ{<!>{p{x)))dx = -E 



In view of p.3p . (|3.2p . in the case where g{k) = l{fc > 1}, the entropy function 
becomes 



Scp_,v+{E) - E = - mi I p{x)\og 



^fA^ _ zmpix))) 

Hpix)) Z{<i>{pix))) ^ ' 

where the infimum is carried over all density profiles p : [0, 1] — > M+ such that 
T^ip-,ip+ip) = Therefore, in the case g{k) = l{k > 1}, where an explicit 

formula is available, up to a linear term, S^__^^{E) is obtained by minimizing the 
free energy functional V^_^ip_^^ over all density profiles p with energy 'Hip_^ip^{p) 
equal to — E'. 

Finally, if we define S^_.tp^ : M — > R by 

SiE) = -inf{]|Vlog|||-log||Mdx:^'[p-p11og$(p-)dx = -i?}, 
we obtain that 

/•i 

S^_.,^4E) = E + + {plog^p)- log Zmp))}dx) . 
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Note that S{E) < and S{0) — 0. As above, we stress that the identity p.4p and 
all formulas thereafter were derived in the case g{k) = l{k > 1}. 

4. The nonequilibrium pressure 

We prove in this section Lemma 12.11 Recall the definition of the function P 
introduced in (j2.6|) . Wc first prove that P is strictly concave. A long and tedious 
computation concluded with the change of variables t = /[x^ + (1 — x)^] shows 
that for 61 7^ 0, 



I A 



If f t \2 



(4.1) 



where A = a^l{oP + (1 - B = /SVi/?" + (1 - /3)«], 



7^ 



and Z(&) is a normalizing constant which makes /ig a probability measure on [A, B\ 
By Schwarz inequality, the first line of the expression of P" without is posi- 
tive. Therefore, P is strictly concave on the interval (— cxd,0). The strict concavity 
on the interval (0, oo) follows from the claim that for all 6* > 0, < a < /3 < 1, 

s{tffig{dt) - ( y" s{t) fie{dt)y < e J t(l-i) (logY^)%e(dt) , 

(4.2) 

It is enough to prove that 

{sit)~ sir)f flg{dr)^lg{dt) < 9 / t(l-t)(log- -) neidt) . 

A<r<t<B J A V 1 — t/ 

Let H{u) = 6'-Hm-^/^ + (1 - uy^/"}. Assume that 6* 7^ 1 and denote by R the 
primitive of given by i?(M) ~ (0 — l)^^{it^~^/^ — (1— u)^~^/^}. Hence, by Schwarz 
inequality, the left hand side of the previous inequality is bounded above by 

t 



{s'{u)f H{u)-^du) ( / H{v)dv) fie{dr)^e{dt) 



A<r<t<l 
B 

du{s'{u)f H{u)~^ I dr I dtmg[r)me{t){R{t) - R{r)} 



1 

Z2 



A 



Since s'{r) = log[r/(l — r)], to conclude the proof of Claim (|4.2p for 9 ^ 1, it 
remains to show that 



A 



dr / dtmg{r)me{t) {R{t) - R{r)} < 1 (4.3) 



for all 6* > 0, 6* 7^ 1 and < A < M < B < 1. 

The left hand side of the previous inequality can be written as 

I^B 2^ Y 

/ rag{t) R{t) dt — / rag{r) dr — / mQ{r) R{r) dr — / mg{t) dt . 

./„ Z J A J A Z ,],, 
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We need to show that this expression is strictly bounded above by 1 for A < u < B. 
Let K he a primitive of mo R and rewrite the previous expression as 

Jas(w) [KiB)-K{u)]M{u) - [K (u) - K (A)] [1 - AI (u)] 
= K{B) M{u) + K{A) [I - M{u)] - K{u) , 

where M(u) = Z^^ J^^ ^■^■mg{r) dr. This expression represents the difference be- 
tween the convex combination of K{A) and K{B), with weights M{u), 1 — M{u), 
and K{u). For A < u < B, this difference is clearly absolutely bounded by the 
variation of K on the interval [A, B] : 

sup |Jy!is(u)| < sup K{v) ~ inf K{v) . 

A<u<B A<v<B A<v<B 

Maximizing over 0<A<i?<l, we get that 

sup |^A.B(it)| ^ sup K{v) — inf K{v) . 

0<A<u<B<l 0<v<l 0<v<l 

A simple computation shows that 

m = -^iog{tv^ + (i-t)V'^}. 

In particular, for 6^1, K{0) = K{1) = 0, K is symmetric around 1/2, K'{\/2) = 0, 
K decreases on the interval [0,1/2] and increases on the interval [1/2,1]. Its 
total variation on the interval [0,1] is therefore equal to —K{l/2) = [9/ [6 — 
l)]log2i-i/^ = log2 < 1. This proves gS]) and therefore Claim g^]) for 9^1. 
The proof for 6* = 1 is identical, the only difference being the explicit expression for 
the primitives. 

The behavior of P in a neighborhood of is obtained through a simple Taylor 
expansion of the integrand. We have 



P{e) = -log2 + 9logh ' 9^-r + 0{9'') , 

Ji 



where 



h = 



dx 



P-aJa ^/x{l - x) 



1 f^[\ogxf + [\og{l~x)]\ 
h = ^TT--. T / , . dx 



and 

" 2(/3-a)7, ^x{\-x) 
This completes the proof of the strict concavity of P on M. 
We now turn to the claim that 

lim {P{9) - 9E^\ = . 
e-i.±oo ^ ^ ' > 

We consider the limit 9 ^ oo, the other one being similar. By definition of P, we 
have to prove that 

lim 9{ log f ^ii^") - E\ = . (4.4) 

e-i-oo \p — a) ' 

A preliminary computation shows that limg^oo log[^(0)/(/3 — a)] = 

The proof of (|4.4p depends on the positions of a and /3 with respect to 1/2, the 
most difficult case being when 0<a<l/2</3<l. Write A(G) as 



^exp{--log[l+(^) J}+ / -exp{-^log[l+(-^) J} 



dx r 1, 

/2 
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We concentrate on the first integral. Since \e'^ — 1 — q\ < (g^/2)el'l, g e M, the first 
integral is equal to 

"1/2 dx 1 /-i/^ dx r ^ X \Si 1 



1 — a; 0.1^ 1 — X 



1-x 



where £{9) is a remainder absolutely bounded by [log2]2 for 6 > 1. The second 
integral in this expression vanishes as oo by the dominated convergence theorem. 
Hence, log[A(6')/(/3 - a)] = + o{0-^), where 0o{9-^) vanishes as 6* t oo. This 
proves (|4.4p . 

We finally consider the last statement of the lemma. By p.7p and by the change 
of variables x ~ , P{0) — 0P'{0) is equal to 

1 riP/i^-f}))" I xlogx- (l + x)log(l + a;) , 

■ dx 



- ^ f'^'^'^ '^'^ u'^logM^ - (1+M^)l0g(l+U'^) 

Wc examine the case 9 f oo, < a < 1/2 < (3, the other ones being simpler. Since 
A{0) converges to a constant as f oo, only the integral has to be estimated. By 
the dominated convergence theorem. 

du u^logw''-(l + u'')log(l + M^) _ 

because the numerator vanishes as ^ oo. On the other hand, the integral in the 
interval — (3)] can be written as 

^/i--^ du log(l + -g-'^) + M-'^ log(l + -g^) 
u{l + u) + 

By the dominated convergence theorem, this expression vanishes as 6* t oo- This 
concludes the proof of the lemma. 

5. Energy band 

In this section, we determine the energy band [E^, £'+], i.e. the range of + §. 
Recall that M is the set of profiles m : [0, 1] [0, 1] and let M+, M~ be the set of 
profiles m of the form m{x) = l{[0,xo]}, rn{x) = l{[xo, 1]}, respectively, for some 
xo e [0,1]. 

Lemma 5.1. For every < a < (3 < 1, 

sup {F(m)+S(m)} = sup {V{m)+S{m)} = E+ , (5.1) 

and the supremum is attained for a unique profile m G 
Proof. By [3] (2.11)], 

V{m)+§{m) = supg{m,F) (5.2) 



where 

r F'(x) -1 

g{m,F) = m(.T) log_F(a;) + [1 — m(.T )] log[l - F{x)] - log } dx 

Jo I 13 -ai 
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and T is the set of all increasing functions F : [0, 1] — > [0, 1] with boundary 
conditions F{0) = a, F{1) = /3. Moreover, the supremum is achieved at the unique 
solution F of the boundary value problem (|2.4[) . 
We claim that for all F in 

sup Q{m,F) ~ sup Q(m,F). (5-3) 

The first identity in (|5.1|) follows from this assertion. To check (|5.3p . fix F in 
J-" and note that the supremum is achieved by m = \{^^xfW S where 
xf = sup{a; G [0,1]; F(a;) < 1/2}, because is increasing. Of course xf = 1 if 
P < 1/2 and = if a > 1/2. 

We claim that the solutions of the variational problem (j5.1|) belong to A^ + . As- 
sume that there exists toq € such that V^(mo)+S(TOo) = sup„jg^{y(m)-|-§(m)}. 
Let Fq be the unique solution of the boundary value problem (|2.4p associated to 
mo, let xq = sup{x G [0,l];Fo(x) < 1/2}, let rrii ~ l{[0,a;o]}, and let Fi be the 
solution of (|2.4p associated to mi. By (|5.2p . 

e(mi,Fo) < a(mi,Fi) = V^(mi) + §(mi) . 

By definition of mo and by (|5.2|) . the previous expression is bounded above by 



sup {y(m) + §(m)} = ^7(mo,i^o) < sup g{m,Fo) = g{mi,Fo), 

where the last identity follows from the argument presented in the previous para- 
graph. Since the first and the last terms in this sequence of inequalities are the same, 
all terms are equal and C/(mi,Fo) = Q{mi,Fi). Since, by (|5.2p . sup^t/(mi,F) is 
uniquely attained at F = Fi, Fq ~ Fi. Therefore, in view of (j2.4p . mi = mo a.s. 

Let TTixg = l{[0,a;o]}, xq £ [0,1], and let F be the solution of the nonlinear 
boundary value problem (|2.4p with m = rUxo- On the interval [0, xq], F{x) = ae"^ 
for some a > 0, and on the interval [xo,l], F{x) = 1 — (1 — /3)e^^^~^-' for some 
A > 0. Since F must belong to C-^dO, 1]), a and A satisfy 

iaae"'-"'' = (1 - /3)Ae'4(^"'"'' , 



1 - (1 -^)e^(i-^o) 



Thus, 



and therefore. 



A 1 



e^(i--«' = ^T-^ , (5.4) 



a-^ Aa a + Al-/3 



-log(A-| + t1°4^A| = 1- (5-5) 

a Va + Aai A Va + Al-l3\ ^ ' 

Moreover, since xq belongs to [0,1], a, A must satisfy max{l — /?, 1 — ae"} < 
a/{a + A)< min{l - a, (1 - (i)e^}. Let /, g : [0, 1] ^ R be given by f{x) = ae^^S 
.g(a:;) = 1 — (1 — (i)e^^^^'''^ and let h{x) = /(x) — 5(2;). Since h is convex and 

< xo < 1, h{xo) = h'{xo) = 0, h{0) > 0, h{l) > 0. Hence, 1 - a < (1 - /3)e^, 

1 - ae" < 1 - /? so that 

1 -/? < — < 1 - a • (5.6) 
a+A- ^ ' 
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By the explicit expression of F and by ()5.4p . (|5.5p . S['mxa) + V['mxo) is equal 

to 

/ log (1^1 da; - (\ogF{x)dx- [ \og(l-F(x))dx 
Jo ^P~ai Jo Jx„ 

= -log (——7-) log « + T^ogf— TTT--?) l°g^ - log(/3-a). 
a \a + Aa/ A \a + Al — p/ 

Therefore, by the concavity of the log function and by (|5.4p . 

+ n-.„) < log{log(-^i) + log(;^Y^)} - log(/3-«) 

By (j5.6p . the previous expression is bounded by 

sup log I log l-^ — ^1 1 - log(/3 -a) = E+ , 



pe[i-^;,i-a] 



where the last identity follows by a direct computation. The last supremum is 
realized for p = 1 — /? if /3 < 1/2, for p = 1/2 if a < 1/2 < /3 and for p = 1 — a if 
1/2 < a. 

Up to this point, we proved that sup,^g^{S(m) + V{m)} < E+. Assume that 
/3 < 1/2 and set xq = 1. In this case, by (|5.4p . a/{a + A) ~ 1 — /3 and all inequalities 
in the previous argument are in fact identities. In particular, §(toi) + V{mi) ~ 
sup^g^{§(m)+V^(m)} = £'+. Moreover, since log is a strictly concave function and 
since by ((5^ a/{a + A) > (1-/3) for < 1, S{mi) + V{mi) > S{mxo) + V{mxo) for 
xq < 1. In the same, way, if a > 1/2, §(mo) + V{mo) = sup„jg^{S(m) + V{m)} = 
E+ and §(too) + V{mo) > S{mxo) + V{mxo) for xq > 0. 

Finally, if a < 1/2 < /3, let x = log(2a)/ log[4a(l - f3)] and observe that 

S(m^) + V{mx) = sup {§(m) + V{m)} = E+ 

and that S(mj,) + V{mx) > S{mxo) + V{mxo) for xq ^ x. □ 
Lemma 5.2. For every < a < /3 < 1, 

inf {y(m) + S(m)| = inf {y(m) + S(m)} = F_ . 

Proof. Recall (|5.2p . Since for each F in J^, F) is a continuous function for the 
weak topology, + S is lower semicontinuous. In view of the explicit expression of 
Q and by Jensen's inequality, y + § is bounded below. Hence, there exists mg G A4 
such that V{mo) + §(^0) — infmGAl{^('7?) + §(m)}. 

Let Fl = inf„ig^{y(m) + §(m)} = inf„jg7vi sup^^^^jr 5(to, F). This expression 
is bounded below by sup^gjr inf^g^ Gim, F). By the explicit expression of G, 

inf Gim,F) = g(?7iF,F) , 

where mp belongs to A^_. If a > 1/2, (3 < 1/2, mp = 1. mp = 0, respectively. 
Otherwise, mp = l{[xp, 1]}, where xp is the unique point where F is equal to 1/2. 

Assume that a > 1/2. In this case, F* > sup^^g^ 0(1, F) = ^(1) +S(1) > Fl. 
To conclude the proof of the lemma it remains to compute sup^gjr ^(1, F) which 
can be done as in the previous lemma. The case /3 < 1/2 is similar. 
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Assume that a < 1/2 < j3. In this case we have that EZ_ > sup^gjr ^/(jtif, F) = 
suPo<^<i supj^g^^ G{l{[x, 1]}, F), where ^ {F e F : F{x) = 1/2}. Fix < a; < 
1. In each interval [0,a;], [x, 1], the variational problem supp^jr^Q{l{[x,l]}, F) is 
similar to the one in (j5.2p . In the interval [0, x] the solution F^fi of this variational 
problem solves the differential equation (|2.4p with m = and boundary conditions 
F{0) = a, F{x) = 1/2. Analogously, in the interval [a;, 1] the solution F^.i solves 
the differential equation ()2.4p with m = 1 and boundary conditions F{x) — 1/2, 
F(l) = (3. These solutions can be computed explicitly and one obtains that 

F.,o(.) = l-(l-a)cxp{-i^^5M^,}, 0<y<x, 

= /3exp{i^^^(2/-l)} , x<y<l, 

sup g{l{[x,l]},F) - a;logi^^^^l^ + (l-x)logi^ - log(/3-a). 
i?e;r^ x l-x 

Maximizing over x we deduce that x* — log[2(l — a)]/ log[4(l — a)/3] is the optimal 
value of X and that 

dx 



E* > loj 



maxjx, 1 — x} ■ 

Moreover, a simple computation shows that for i^^* o (*''*) ~ ^x* Hence, 
the function G : [0,1] [a, (3] defined by G{y) = F^,,o{y)My e [0,a;*]} + 
Fx* ^i{y)l{y e {x* , 1]} belongs to F and solves the boundary value problem (|2.4p 
for m = l{[a;*, 1]}. Hence, by definition of and by (|5.2p 

b; < S(l{[xM])) + = 5(l{|x-,ll),G) 

^■o.{-^r ,\ , }. 

L p — a Jq, max|a;, 1 — x} J 
This proves the lemma and shows that a profile with minimum energy is given by 



log[2(l - a)] ^11 
log[4(l - a)/3] ' JJ- 



□ 



We proved in the previous lemma that 



inf sup Q{m,F) ~ sup inf Q{m,F) . 

Fix < a < /? < 1 and assume that V{m) + S{m) = Then, 

m(l - m) = a. s. (5.7) 

Indeed, fix a profile in Ai. Since 1/(to) + §(to) = Q{m,F), where _F is the 
solution of (|2^ . 

J[y + §](m) _ (5^(771, F) 5a(m,F) (5F 
5m 6m 5F 6m 

Since F solves the Euler equation (|2.4I) associated to the variational problem 
supG^;(77i,G), 5G{m,F)/6F = 0. On the other hand, 50/ 5m = -log[F/(l - F)], 
so that 

5[y + S](m) , F 

— ^ — — — - = — log • 

5m ^1-F 
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This expression does not vanish because F is strictly increasing. Therefore, the 
extremal values of y + S are attained at the boundary. 

This formal argument can be made rigorous. By the proof of Theorem 7.1 in 
[6] , V and therefore 1/ + § is Gateaux differentiable, and the Gateaux derivative of 

+ § at m is equal to - log[F/(l - F)]. 

It follows from the previous results and the variational formula p.Sp that for 
< a < ^ < 1, 

S{E) = -oo for E ^ [E^,E+] and S{E±) = 0. (5.8) 

6. ISENTROPIC SURFACE 

We determine in this section the isentropic surfaces defined by 

£k = {E>0: So^AE) = A-}, K G [0,log2] . 

6.1. The equilibrium case. Assume that a = f3. We have already seen right 
after (|2.3p that the energy band is reduced to the point log 2 in the case a = 1/2. 
Assume therefore that a ^ 1/2 and fix K e [0,log2). There exist exactly two 
solutions < m^{K) < 1/2 < m+{K) < 1 of -s(m) = K. Hence, in view of p.3p . 
the level set 8k = {E^{K), E+{K)}, where 

E±{K) = - log(l - a) - log- m±{K) . 

I — a 

For K = log2, £k is the singleton { — (l/2)[loga + log(l — a)]}. 

6.2. The nonequilibrium case. Assume now that a < (3. Let D : R — !• M be 

given by D{0) = 0P'{e) - P{0). Since D'{0) = eP"{0) and since, by Lemma O 
P is strictly concave, D is strictly increasing on (— oo; 0] and strictly decreasing 
on [0, +oo). Moreover D{Q) = log 2 and, by Lemma [2.11 \\m0^±Qc, D{9) = 0. In 
particular, for every K G (0,log2) there exist exactly two values 9-{K) < < 
e+{K) such that 0±iK)P'{e±{K)) - P{0±{K)) = K. 

Fix < /\ < log 2. By ((2T3)) . E belongs to £k if and only if K = 0eP'{Oe) ~ 
P{0e) = D{eE), where P'{0e) = E. Hence Oe = 0±{K) so that E = P'{e±{K)) 
and 

£K = {P\0+{K)),P'{e^{K))] . 

If we let a and /3 vary, we see that the i^T-isentropic surface is composed of the 
two manifolds E^ = P'{0+{K)) and Ej^ = P'{0^{K)) which satisfy 

E-^{a,P) < E+{a,(3) . 

7. Comparison with local equilibria 

In this section, we compare the entropy function Saj3 with the entropy function 
associated to product measures with a slowly varying density profile that will be 
called local equilibrium entropies. 

Let p be the product probability measure on {0, 1}^^^ given by 

N-l 

<piv) - l[p{^./Nr^{l-p{x/N))'~^^ , 

x=l 
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where p : [0,1] is the stationary profile p{x) = (1 — x)a + xp. Denote 

by S := Sa.p ■ R+ — > {^oo} U [0, log2] the entropy function corresponding to the 
Gibbs state p'- 

S{E) = hm hm llog ^ 1 {\N-Hog,^^^f,{rj) + E\ < 6} (7.1) 

whenever the hmits exist. 

Let P := Pa.i3 : R — M be the function defined by 



1 /■'^ / 1 



Lemma 7.1. P is a strictly concave function and 



dx 



lim ^ = E^, hm - eP'iO)} = 0, 

w/iere 

1 1 ~ 1 1 

^- " « — / — r~i — T ' ^+ = « — / ■ r 1 — X dx . 

p — aj^ niax|a;, 1 — xj p — aj^ min|2;, 1 — x} 
Moreover, as 9 0, 

p[e) = -iog2 + /^og ^ — =dx + o{e') . 

The proof is elementary and left to the reader. It follows from this result that 



lim P\e) = E^ 

I— >±oo 



Proposition 7.2. For E >0, 



S(E) = inf \0E-P(e)] . 

If E belongs to the energy band (£'_,i?+), the infimum is attained atOs = 9e{ol^I3) 
the unique solution of P'{9) ~ E and S{E) = OeE — P{9e)- S{E) = — oo if 
E i {E^,E^\ and S{E±) = 0. 

Proof. Multiplying and dividing the indicator in (j7.ip by 2^~^, we reduce the 
computation of the entropy to a large deviations problem for independent Bernoulli 
random variables and we obtain that 

S{E) = sup {§(m) : /(m) + §(m) = E} , 

rriGM 

where / stands for the large deviations rate function given by 

/"^ r / \ 1 ^{x) 1 — 771(3::) ~i 

lim) = / < mix) log _. , + 1 - m(x)\ log ——- > dx . 

Jo P[x) 1 - p{x) i 

One should compare this expression with the variational formula (|2.5p for the 
nonequilibrium entropy. 

Repeating the arguments presented in the proof of (|2.9p . we deduce that 
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where 9 is the unique solution of P'{d) ~ E. The rest of the proof is similar to the 
one of ([21^ . □ 

Let So = ^"(0), ^0 = ^'(0). By Lemma O and O 

Eo = log (-^ , . dx) , ^0 = f log ( , .| J rfx . 

By Jensen's inequality, -E- < i?+ < and Eq < Eq- Since min{x, 1 — < 
1/2 < max{a;, 1 — a;} and y^x{l — x) < 1/2, we may compare all variables with log 2 
to obtain in the end that £'_ < £'_ < log 2 < Eq < minji^o, E+} < max{i?o, E+} < 
E+ in the case a < (3. 

The nonequilibrium and the local equilibrium entropy differ. For every < a < 
^ < 1, 5 < 5 in the interval {E-,Eq) and 5 < 5 in the interval (i?o, E+). 

Indeed, fix e {E^,Eq). By Jensen's inequality, e-^P{0) < 0-^P{e), 6* e M. 
Therefore, for every 9 > 0, 9E - P{9) < 9E - P{9). On the other hand, since 
E < Eq and 9eo = 0, P'i9E) = E < Eq ^ P'i^Eo) = P'{0)- Hence, since P' 
decreases, 9e > 0- A similar argument shows that 9e > 0. In conclusion, by the 
variational formula for the entropies presented in (|2.13p and [721 

S(E) = M{9E~P(9)] = M {9E ~ P(9)} 

< M{9E-P(9)] = M{9E-P(9)] = S(E) . 
A similar argument proves the other claim. 

Consider the sequence of random variables Yiq{vi) = —N~^ log ^^^^{r]), Y^ir]) = 
—N^^ log i^apiv) defined on the probability space fi^v equipped with the probability 
measure fJ.^ i^afi^ respectively. By (|2.ip . the sequence {Yjy : N > 1) satisfies a 
large deviations principle with convex rate function J{E) = E — S{E). By similar 
reasons, the sequence (Yn '■ N > I) satisfies a large deviations principle with convex 
rate function J{E) ^ E- S{E). 

Bahadoran [T] proved that E^n ^[Y^] and E^n ^[Y^^] have the same limit given 
by the Gibbs-Shannon entropy 



E ^ — s{p{x)) dx , 
Jo 

where p{x) = a + (/3 — a)x. This result can be recovered from ours. 

Lemma 7.3. The nonnegative rate functions J and J vanish at the same and 
unique point 

E ^ P' {1) = P' {!) . 
In particular, Y^r under p, and Yn under p converge in probability to E. 

Proof. By Lemma [2.11 the variational formula (|2.13p and the assertions following 
this formula, the nonnegative rate function J is strictly convex on [E^,E^], dif- 
fercntiable in {E-,E+), and equal to +oo outside of the interval It has 

therefore a unique minimum E G [£'_,i?+]. 

By ([2J3l) and ([2S|), S'{E) = 9e on {E^,E+), where 9e is the unique solution 
of P'{9) — E. By Lemma [2?T| limg^ioo -P'(^) — Ez^z. It follows from the previous 
two facts that \uiie^e± J'{E) = ±oo. Since J is strictly convex, J has a unique 
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minimizer E in {E-,E+) solution of 6'^ = S'{E) = 1. Applying P' on both sides 
of this equation, we deduce that E = -P'(l)- 

We claim that J{E) = 0. To prove this identity we need to show that S{E) = E 
or, in view of (|2.13p . that 9j^E — P{6j^) = E. Since = 1, this equation is reduced 
to P(l) = 0, which is easy to check in view of the explicit formula (|2.6I) for the 
noncquilibrium pressure. 

The same argument applies to J and the result follows from the identity 



P'(l) = P'(l) 



s{a + (/3 — a)x) dx 



□ 



In [SJ Section 7] , the authors compute the limit of the variance of the sequences 
(Yat : iV > 1) and (Y/v '■ N >1) and show that the limits differ. This result can be 
recovered from a second order expansion of the entropy function Sa,fi. 

Wc have seen that the rate function J has a unique minimum a,t E. It is well 
known from the theory of large deviations that the asymptotic variance of the 
sequence Y/v is given by J"{E)~^ = —S"{E)~^. Since = 1 and since S is 
the Legendre transform of the nonequilibrium pressure P, we have that S"{E) = 
1/P"{9j^) = 1/P"(1). Hence, — P"(l) is the asymptotic variance of the sequence 
Yn. 

By taking 6* = 1 in ()4.ip we obtain that 

-P"(l) = /%(l-t)[log 



/3~a 
1 



2(/3-a)2 



1 -t 



dt 



/3 



dx 



13 



(7.2) 



dy{s{x) ~ s{y)f 



A long and tedious computation shows that this expression coincides with the 
limiting variance derived in [9]. 

A similar computation in the equilibrium model gives that the asymptotic vari- 
ance of the sequence Yat is equal to 



t{l-t) 



loa 



1 - t 



dt 



In particular, the asymptotic variance in the nonequilibrium model is strictly bounded 
above by the asymptotic variance in the equilibrium model. 



Appendix A. The assumption (H) 

The stationary state /^^^ of the symmetric simple exclusion process with open 
boundary conditions can be expressed in terms of a product of matrices [9] : There 
exists matrices D, E and vectors \V), {W\ such that 

DE - ED = D + E , {{l-f3)D- f3E} \V) = \ V) , 
{W\{aE-{l-a)D} = {W \ 

and 



{W\{D + E)N-i\V) ' 



(A.l) 
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where the weight uiNiv) is given by 



N-1 



^Niv) - {w\ n M^)D + [1 - vi^m \V) . 

x = l 

The partition function {W\{D + E)^-^\V) can be computed exphcitly [H (3.11)]: 
{W\iD + E)''-'\V) = ^f^^^- (A.2) 

Decompose the set {1, . . . , — 1} into r — adjacent intervals Ki, . . . , Kr 
of size eN and denote by M = {Mi, . . . , Mr) the number of particles in each box. 
We recall that • |M) is the probability measure /i^^ conditioned to have Mj 

particles in Kj, j = 1, ... ,r. 

Let r]^ , rj~ be the configuration in fiAr(M) with all particles in each interval Kj 
at the left most, right most positions, respectively. Hence, if Kj = {xi, . . . 
Mj — M, rf^{x) = 1 if and only if < a; < xm, rj^ {x) = 1 if and only if 
XL-M+i < X < XL- 
Lemma A.l. For < a < P < 1, ij e riAr(M), 

fi{r]^\M) < fi{T]\M) < fi{i-i+\M) . 

Proof. This is a simple consequence of the matrix product form of the stationary 
state. Let 77 be any configuration in fljv and let 1 < a; < — 1 be any site such 
that 77(2;) = 1, rj{x + 1) = 0. Then, uj^io-'^'^^^'ri) ~ ^n{v) is equal to 

an-l Af-1 

{W\l[My)D+[l-r^{y)]E}[ED-DE] {ri{y)D + [1 - v{y)]E} \V) 

y=l y=x+2 

x-1 N-1 

= -{W\l[My)D + [l-rj{y)]E}[D + E] J] {ri{y)D + [1 - 7M]E} \V) 

y=l y=x+2 
= -UJN-liO - '^N-liO < 0, 

where ^, C are the configuration of flN-i given by ^ = (f?(l), • . . , r]{x — 1), 1, r]{x + 
2), . . . , r,{N - 1)), C = Ml), ...,vix- 1), 0, r,{x + 2),..., r,(7V - 1)). □ 

Hence, the derivation of the assumption (H) is reduced to the proof that 
1 



lim lim sup sup 

e^O N^oo M N 



logM^:,^(r/+|M)-logM^:,^(r7-|M) = 
or, equivalently, to the proof that 

' 







lim lim sup sup — 

For each fixed e > and M, the configurations r]~^, rj~ are associated to density 
profiles p+, p- : [0, 1] ^ [0, 1] defined by p+ = Ei<,<. " " l)e + Pi)}, 

p^ = Si<i<r -i{[*^ ~ Pi, *^)}' ■where pi = Mi/N, 1 < i <r and 1{A} stands for the 
indicator of the set A. Therefore, by the large deviation principle for the density 
profiles under the stationary state pi [71 [TT|. 

1 

n'-^'oo N 



lim — log^^^^(?7±) = -Va^fiip^) , 
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where Va,p is the functional introduced just before ()2.4p . 
Hence, to prove assumption (H), it remains to show that 



hm sup 



= 



where the supremum is carried over all < pi < e, 1 < i < r. Since p^{x) is either 
or 1, we may replace in the previous formula, Va.p by Va,i3 + §■ By [3l Theorem 
4.1], this functional is continuous in A4. For each £ > 0, denote by p^'"^ the profiles 
which attain the previous supremum with Va,i3 by Va^p +S. By compactness of A4, 
there exists a subsequence Sk iO for which 



lim sup 



v^Ap^-'') + Hp^^"') - v^Ap-^'-) - s{p--'-) 



= lim 

and converges weakly to some profile p. Clearly, the sequence p"'^'' converges 
weakly to the same profile p. Since Va,i3 + S is continuous in Ai, assumption (H) 
is proved. 
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